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ABSTRACT 

\I'or  continuous  observations  fi'om  time-sequential  studies,  suitable 
Cramer-von  Mises  and  Kolmogorov-Smirnov  type  (nonparametric)  statistics 
(based  on  linear  rank  statistics)  for  testing  hypotheses  on  some  multi- 
ple regression  models  are  proposed  and  studied.  Asymptotic  theory  of 
these  tests  is  provided  for  both  the  null  and  (local)  alternative  hypo- 
theses situations  and  is  based  on  the  weak  convergence  of  suitable  rank 


order  processes  (on  the  0 [ 0 , 1 ] 
ian  motions.  Bahadur  efficiency 
values  of  the  percentile  points 
posed  test  statistics,  obtained 
provided. 


space)  to  certain  functions  of  Brown- 
results  are  also  presented.  Empirical 
of  the  null  distributions  of  the  pro- 
through  simulation  studies,  are  also 
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1.  INTRODUCTION 


In  longitudinal  (time-sequential)  studies  relating  to  clincial 
trials  and  life  testing  problems,  the  experimenter  sets  out  to  plan 
beforehand  as  to  the  maximum  number  of  responses  to  be  observed  or  the 
maximum  duration  of  the  experimentation.  As  the  experiment  is  continu- 
ously observed  over  time,  even  with  this  restricted  design,  the  experi- 
menter has  the  option  of  reviewing  the  outcome  as  the  experiment  pro- 
gresses, enabling  him  to  terminate  the  experiment  at  an  intermediate 
stage  if  the  cumulative  evidence  indicated  by  the  data  at  that  stage  is 
strong  enough  to  reject  the  null  hypothesis  and  further  continuation  of 
experimentation  can  nut  lead  to  a different  inference.  This  pseudo-  se- 
quential test  procedure  (which  is  distinguished  from  the  classical  sequen- 
tial test)  arises  from  what  is  called  progressive  censoring  schems  (PCS) 
as  at  the  successive  censoring  time-points  (responses  or  failures),  the 
test  statistics  are  based  on  uncensored  subjects  only.  It  may  be  men- 
tioned that,  for  the  applications  of  this  procedure,  it  is  not  necessary 
that  one  works  with  a restricted  design  only. 

For  two-sample  location  and  scale  problems  as  well  as  the  simple 
regression  model,  Chatterjec  and  Sen  (1973)  have  developed  a general 
class  of  (linear)  rank  statistics  incorporated  for  testing  under  PCS. 

In  the  current  investigation,  their  theory  is  extended  to  the  multiple 
regression  problem  which  includes  the  multi-sample  location  problem  as 


a special  case, 
here. 


Also,  a wider  class  of  test 
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Section  2 is  devoted  to  the  basic  formulation  of  the  problem.  Sec- 
tion 3 deals  with  the  development  of  PCS  tests:  two  different  types,  viz., 
Kolomogorov-Smirnov  and  Cramer-von  Miscs,  both  based  on  linear  rank  sta- 
tistics, are  considered.  Sections  4 and  5 are  devoted  to  the  study  of 
the  asymptotic  distribution  theory  of  the  test  statistics  under  the  null 
and  (local)  contiguous  alternative  hypotheses.  Asymptotic  (Bahadur-) 
relative  efficiency  results  are  presented  in  Section  6.  Section  7 deals 
with  the  simulation  study  of  the  null  distributions  of  the  proposed  test 
statistics.  In  Section  8,  comparison  of  the  proposed  PCS  tests  with  the 
fixed-plan  censoring  modifications  of  the  Kruskal-Kallis  test  by  Basu 
(1967)  has  been  made  and  the  scope  of  applicability  of  the  PCS  tests  in 
the  context  of  right  truncation  is  discussed. 

2.  PRELIMINARY  NOTION'S 

Let  (X^,  iei}  be  a sequence  of  independent  random  variables 
with  continuous  distribution  functions  (d.f.)  {F^,  i>l},  specified  by 

the  model 

F.(x)  = F(x-30-6'c.)  , - oo  < x < 00  , is]  , (2.1) 

where  in  this  conventional  multiple  regression  model,  the  d.f.  F is  not 
known,  3 , B'  = (B. , . . • , B ) are  unknown  parameters,  p(s  1)  and  (c.  ' = 

U 1 P "^1 

(Cj.,...,c  isi}  is  a sequence  of  (known)  vectors  of  regression  con- 
stants. Our  concern  is  to  test 

Hq:  3=0  vs.  Hj : fi*  0 , (2.2) 


L. 
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(trcnting  3^  as  a nuisance  parameter  and  without  assuming  the  form  of 
F to  be  specified),  when,  in  fact,  we  have  a life  testing  model,  as  may 
be  posed  below. 

For  every  N(>  1),  let  Z^j  $ ...  i be  the  ordered  random  vari- 

(N)  ’ 

able?  corresponding  to  (X^,...,X^)  =X  , say;  by  virtue  of  the  assumed 
continuity  of  the  F^,  tics  among  the  observations  may  be  neglected,  in 
probability.  Let  SN  = (RN1 . • • • ,RMV)'  and  £w  = (SM1  , • • . ,SMM)'  be  respec- 


NN 


-N  NI : 


NN 


(N) 


tivcly  the  vectors  of  ranks  and  anti-ranks  of  the  elements  of  Xv  ' , so 
that  (tics  neglected) , R.,g  - S^.R  = i,  l<i<N,  X.=Z^_  , lsiSN.  In 

a life  testing  problem,  one  typically  observes  the  successive  order  sta- 
tistics {Z^}  along  with  the  corresponding  {S^},  and  based  on  a part 
of  the  sequence  {7  , S^;  IsisN},  the  problem  is  to  test  for  Hq  in 
(2.2).  The  PCS  test  is  a pseudo-sequential  procedure  where  a early  termi- 
nation of  experimentation  is  feasible  if  observing  (Z.,.,  S...,  1 ^ i s k} 
for  some  k(s  N) , the  accumulated  statistical  evidence  leads  to  a deci- 
sive conclusion.  Our  proposed  FCS  tests  are  based  on  suitable  linear  rank 
statistics  which  we  introduce  in  this  section.  The  actual  test  statistics 


will  be  introduced  in  the  next  section. 

For  every  N(2  1),  wc  conceive  of  a set  of  scores  (a^(l) , . . . ,a^(N)) 
generated  by  a score  function  <t>  = (<{)(u),  uc  (0,1)}  in  the  following  way: 
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Further,  4>(u)  is  assumed  to  be  expressible  as  the  difference  of  two 
[(J)j(u)  and  (^(u)]  non-decreasing  and  square  integrablc  (inside  I - 
[0,1])  functions.  Note  that  (2.4)  holds,  in  particular,  when 


aj,(i)  = E4»(UNi)  or  «J»(i/(N+l))  , l<i<N, 


(2.5) 


where  ll^j  < ...  < are  the  ordered  random  variables  of  a sample  of  size 
N from  the  rectangular  [0,1]  d.f.  In  order  to  simplify  the  notations, 
in  the  sequal,  we  let 


= 1' 


(p  = I <t(u)du  “ 0 , A 
0 


= jUi-lV1^0  and  A,5=  (N-D-'^it^a)-^]^1,  N>2  . (2.7) 

Now,  at  the  k-th  stage,  the  observable  random  variables  are  {Z^, 

S^.  ; lsisk},  and  based  on  these,  we  define 


♦ 


<[>^(u)du  - = 


-lvN 


2 

<J>  (u)  du  = 1 , 

0 


...  - ,2 


(2.6) 


.(1) 


Vk=  (^ksN) 


(p),, 


N,k' 


by  letting 


^,k  = ^i=l(£S  ‘%PaN^)  + j^i=k+l(~SN."£N^(N'k)  ^i=k+laN(l)^ 
Ni  *-  Ni  -1 


(2.9) 


where  c„  = N ,c.  and 
~N  Li= l~i 

(N-k)’1^=k+1aN(i)  , lskSN-1  , 
[0  , k = 0,  N , 


a*(k)  = \ 


(2.10) 


and,  conventionally,  wc  let  TM  = 0.  Note  that  Tx.  M . =T.,V  =T.f,  where 

~N , 0 ~ ~i\  , N - 1 ~NN 


Tm  = yN  ,(c„  -c.,)av(i)  = yN  (c.  -cM)aw(RM. ) . 

~N  ^l-l  ~SkI.  ~N  N ^1=1  ~i  ~N  N Ni 

Actually,  If  we  let 


Ni 


bN,kW  = 


aN(i)  , 1 * i s k , 

a*(k)  , k < i ^ N , 


k = 0,1... .,N  , 


then,  we  may  rewrite  k as 

^N,k  = ^i=l(~i"%-)bN,kfRNi-)  * k = 0,1,‘,-,N  • 


(2.11) 


(2.12) 


(2.13) 


Note  that  by  (2.7),  (2.10)  and  (2.12),  for  every  l<k<q<N, 
(N-l)-1^=1bNjk(i)bN>q(i)  = (N-ir^jaJtf)*  (N-k)(a*(k)]2} 

" ^-(N-D'^i^CD-a^fk)]2  (2.14) 


= 1 - (N-l)"1^=k+1[aN(i)  - a*(k)]2  = A2>k  , say  , 

where  A2jk  is  in  k (0sk<N)  and  A^Q  = 0,  A^N1  = =A^  = 1 . 

Let  us  also  denote 


CN'^1(£i-cs)(c.-cN)'=  


(2.15) 


and  assume  that  there  exists  a positive  number  N^,  such  that 


is  positive-definite  (p.d.),  V N 5 Ng  . (2.16) 

We  may  remark  that  under  in  (2.2),  = (Z,^ , . . . ,Z^)  anc* 

(and  S^)  are  stochastically  independent;  (and  S^)  assumes 

all  possible  permutations  of  (1,...,N)  with  the  equal  probability  (N!)  *. 
Hence,  from  (2 . 13) - (2 . 15) , we  obtain  by  some  routine  computations  that 


T 
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E(3^,k'V  = ~ ’ 0;Sk-N  » 


E(T  ,T'  |tl  ) = Aki  , • Cv  , 0<k  , q<N  , 

~N,k~N,q'  0 N,kAq  ~N  > m > 

where  aAb  = min(a,b).  Finally,  we  define 

fo  , k = 0 , 


k 


N , k 


[(~N,k~N  ~N,k)  ' k=:  k - • • ,N  . 


Our  proposed  tests  (under  PCS)  are  based  on  the  partial  sequence 


{LN,k:-°-k-r} 


where  r is  any  pre-assigned  integer,  such  that 


(2.17) 

(2.38) 


(2.19) 


(2.20) 


r/N  -*■  6 : 0 < 6 £ 1 , as  . 


(2.21) 


In  reality,  mostly  6 is  less  than  one. 


3.  PCS  RANK  ORDER  TESTS  FOR  NO  REGRESSION 

For  every  N>ri  1,  we  introduce  a sequence  (k^  r(t),  t e 1}  of 
interger-valued,  non-decreasing  and  right-continuous  functions,  where 

kN  r(t)  = max{k:  A^  ^ ~ tAN  ’ t £ 1 * C3-1) 

2 

and  A^  ^ is  defined  by  (2.14).  Consider  then  the  process  k r = ^N  * 
tel}  by  letting 

V(t)  ■ V!rLN,kv  (.)  • <3-2> 

N , r 

Thus,  Ym  belongs  to  the  space  D[0,1],  endowed  with  the  Skorokhod  J.- 

* * i r 1 

topology.  IVe  propose  the  following  two  test  statistics: 
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(i)  Kolmogorov-Smirnov  (KS-)  type  statistics.  Wc  define 


^N.r  = tel^N.r^  ^ 

- / max  t'  C'^T  1 

N,r\o<k<r  ~N , k~N  ~N,kJ  / ‘ 


Note  that,  if  we  let 


C ■ 0£ksr  • 

then,  by  definition, 

0 = < ...  <K^r)  = K* 

W,r  N,r  N,r  N,r 

(ii)  Cramer-von  Mises  (CvM-)  type  statistics.  Here,  we  define 


(3.3) 


(3.4) 


(3.5) 


‘J.r  ■ J0  VLtt)dt 


H*.  = J,  »*l,„ /.c 

N,r  *^=0  k Nk  .V, r 

= a'2  irhtir  .c:1!., .) 

N.r^kO  k ~n,k~N  ~N,k 


(3.6) 


If,  we  let  t^  = A2  k/A2  r>  0<k<r,  and  define 

Xk“tN^1)"tS"AN!rI(N-k-1)/(N-1)(N'k)llaN(k+1)^  > (3-7) 

O^lcfiN-l,  then  from  (3.1),  (3.2),  (3.6)  and  (3.7),  we  have 

_ rr-l,*,2  .,2 


(3.8) 


where  the  A*  also  depend  on  N,r.  Here  also,  if  we  let 

.00 


0 < k < r , (3.9) 


wc  obrain  that 


1 
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The  monotonicity  property  in  (3.5)  or  (3.10)  is  then  incorporated 
in  the  formulation  of  the  following  PCS  tests.  For  a preassigned  level 
of  significance  ot(0<Ci<l),  let  K*  , and  W*  be  defined  by 


,r,cJV  ; 

(3.11) 

(3.12) 

Then,  for  N items  under  life  testing,  as  the  successive  failures  are 

(k1  fkl 

observed,  at  each  failure  Z^j.,  we  compute  ^ (or  Nl^  p,  k ^ 0.  If, 
for  the  first  time,  for  some  k(<  r) , is  > k*  (or  is 

> Mj  ),  experimentation  is  discontinued  following  Z , along  with  the 
rejection  of  11^.  If,  no  such  k(5  r)  exists,  experimentation  is  curtailed 
following  along  with  the  acceptance  of  H^.  Note  that,  by  definition, 

K*  and  M*  are  both  functions  of  R..  (and  c,,...,c..),  and  hence, 
under  in  (2.2),  they  are  distribution-free  statistics.  Thus,  both  the 

PCS  tests  based  on  and  {M^^}  are  genuinely  distribution-free 

(under  H^) . Also,  by  (2.13),  (2.19),  (3.3)  and  (3.8),  these  statistics 
remain  invariant  under  any  non-singular  transformation  on  the  regression 
vectors  {c.}.  That  is,  if  we  let  d.  = c.  + Tc.,  i > 1 , where  c_  is  arbi- 
trary  and  T is  non-singular,  and  for  (2.1),  we  rewrite  x-Bq-B'c^  as 
x-yQ-y'd^,  lsiSN,  where  y'  = B'T  1 and  y0  = 8q-y'cq>  then  repla- 
cing in  (2.11),  (2.13)  and  (2.19),  the  c^  by  d^  and  denoting  the  result- 
ing  statistics  by  L^,  we  have  ^ V k>0  and  £ (non-singu- 

lar). Since,  in  many  cases  (vis.,  the  multipsamplc  location  problem),  the 
c,.  and  B in  (2.1)  are  not  uniquely  defined,  this  invariance  is  rather 


L 
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important  and  it  eliminates  the  arbitrariness  of  the  choice  of  the  c. 

~i . 

Thus,  the  proposed  PCS  tests  are  invariant,  distribution-free  tests,  for 
small  values  of  N,  the  exact  null  distribution  of  ^ r or  r can 
be  derived  by  direct  enumerations  of  the  exact  distribution  of  (or 

S ) over  the  M!  equally  likely  permutations  of  (1,...,N)  giving  rise 

rrl 

to  N equally  likely  realizations  of  (S^, . . . ,S^  ) from  (1,...,N). 
The  task  becomes  prohibitively  laborious  as  r (and  hence,  N)  increases; 
for  this  reason,  we  take  recourse  to  the  asymptotic  distribution  in  the 
next  section. 


4.  ASYMPTOTIC  NULL  DISTRIBUTION  THEORY  ' 


To  study  the  asymptotic  distributions  of  K*  ^ and  M*  ^ (both  of 

which  are  functionals  of  the  process  Y , defined  in  (3. 1) - (3. 2) ) , 

n , r 

first.  We  consider  the  weak  convergence  of  {Y^  to  appropriate  func- 
tionals of  Browr.ian  motions.  Let  W_.  = {W^  (t) , tel},  j = l,...,p  be  p 
independent  copies  of  standard  Brownian  motions  on  I,  and  define 
Y = (Y(t),  tel}  by  letting 


Y(t)  = [lj  = 1W- (t)]  2 . tel. 

Then,  Y belongs  to  the  space  C[0,1]  with  probability  1. 
stage,  we  introduce  the  following  (Noether-type)  condition: 


(4.1) 


At  this 


lim  s up f max  - , -1 

N-»~  l<i<N '■~i  ~NJ  ~N 


(c.-c  ) 
~i  ~N  \ 


0 . 


(4.2) 


Then,  we  have  the  following 
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Theorem  4.1.  Under  11^  in  (2.2),  (2.21),  (4.2)  and  the  assumptiom 
made  on  the  scores  in  Section  2,  as  N + °°, 

Y ^ Y , in  the  J -topology  on  D [0 , 1 ] . 

, T 1 


(4.3) 


Proof.  We  have  noticed  in  the  last  section  that  the  L . (and  hence, 

, K 

Y^  r)  remain  invariant  under  any  non-singular  transformation  on  the  c^ . 
Thus,  under  (2.16),  there  exists  a non-singular  matrix  JF.^  such  that 


E’CnFn  = Ip  = .Diag(l , . . . , 1)  . 


Let  then 


r-.  ~ r-  ^ ~ 

c.  = Lt(c-  '£».)  . i - i - N»  so  that  ) £.  = 0,  Z = i 
~i  ~N  ~i  ’ ’ . Lj-  ~i  . u ~i~i 


i=l 


i=l 


~p 


(4.4) 


(4.5) 


Further,  in  (2.13),  replacing  c.-£N  by  £ i , l<i<N,  the  resulting  vec- 
tor is  denoted  bv  T , for  k = 0,l,...,N.  Then,  by  (2.19),  (4.4)  and 
' S , k 


(4.5),  we  have 


V • for  k*»-' N • 


(4.6) 


Defining  {k  (t),  tel)  by  (3.1),  we  introduce  a p-variate  sto- 
, r 


chastic  process  h’N  r = (WN  r(t),  tel)  = {(W^  r(t) , . . . .^^.(t))  , tel) 
by  letting 

V-w  ■ *i!A,iN,r(t)'  teI-  (4-7) 

Then,  by  (3.2),  (4.6)  and  (4.7),  we  obtain  that 

Yr(t)  " 


t € I 


(4.8) 


Hence,  i f we  let  W*  { (W.  (t) , . . . ,W  (t) ) ' , tel)  where  the  W.  arc 

~ . 1 P J 
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defined  prior  to  (4.1),  then  to  prove  (4.3),  it  suffices  to  show  that 


in  the  J -topology  on  D[0,1]  , 


(4.9) 


and  for  this,  we  need  to  show  that  (a)  the  finite-dimensional  distribu- 
tions (f.d.d.)  of  { IV. . } converge  to  those  of  W and  (b)  (W.,  } is 
~N,r  ~ ~N,r 

tight. 

Note  that  under  the  assumed  regularity  conditions  on  the  scores,  by 
Lemma  4.3  of  Chatterjee  and  Sen  (1973)  and  our  (3.7), 


max 


0<lc<r  k 


X*  0 as  N 00  . 


(4.10) 


(.K  „(t)) 


^ N J 

As  a consequence  , liin^  t^  y = t for  every  t belonging  to  £0,1J. 

Hence,  if  for  q(s  1)  and  Ostj  < ...  <t  i 1,  t = (t^, . . . ,t  ) ' , then 

by  (4.4),  (4.5),  (2.17)  and  (2.18),  we  have  on  letting  W,  (t)  = 

, r 


ttllr'V’- 


’~N,r  q 


(t  )]  that 


%,,.(!) -2  and  V[BN  r(t)]  -IpS  C(t.«tp) 


(4.11) 


where  V[  ] stands  for  the  dispersion  matrix  (of  order  pq * pq)  and  ® 
for  the  Kronecker  product.  Defining  W(t)  = [ W* (t^) , . . . ,W' (t  ) ] , it  fol- 
lows by  routine  steps  that 


EW(t)  = 0 and  V[W(t)]=X  ®((t.At„))  . (4.12) 

~ ~ ~ ~ ~ ~ -wp  J C 

Thus,  to  prove  the  convergence  of  f.d.d.'s  of  (W^  t.o  those  of  W , we 

need  to  show  only  that  for  any  given  q(s  1)  and  t=  (t  ,...,t  )',  W (£) 

~ 1 q in  , r ^ 

is  asymptotically  normal.  Writing 
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Vr®  ■ W * (kl V'  1 (4-'3) 

and  noting  (4.7)  and  (4.10),  it  suffices  to  show  that  (T^  ^ ^ • »XN  j.  ) 

has  asymptotically  a multinomial  distribution.  We  consider  the  rolled 
out  vector 


linear  compound 


;(D 

T (p) 

?(1) 

N, kj  ’ 

•"’N.y 

•’  N,k 

(gjl 

’•••’gpl”" 

,Rlq’ ‘ 

= S’li: 

rN  V] 

) ~ 

M,q  _ H=1 

n=  1 SnS, 

i 

15), 

(4.5)  and  (4 

• 15), 

isiski ; 

. ...  l Vv) 

JN(i)  = • ^ L J v=kj  + l J 


for  ks  ] < i < ks,  s = 2, . . . ,q  ; 


r i N 

I11  1!°'  lS  • I aM(v)  , k <i<N  . 

1 j =1 pmj  N-k.  ^ +1  N q 

L j J 


Thus,  if  we  let 

f o = H ,0-  for  i = 1 , . . . ,N  , 

SNi  m"1  mSNi 

we  have  from  (4 . 15)  - (4 . 17)  that 

u;<s>  ■ • 


where  under  1L  in  (2.2),  = (SN1  > • • • .SNN)'  assumes  all  possible 
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permutations  of  with  equal  probability  (N!)  Thus,  U*(g) 


is  a simple  linear  (anti-)  rank  statistic.  Note  that  by  assumptions  (2.3)- 


(2.7)  and  by  Lemma  2.2  and  Theorem  3.2  of  Hajek  (1961),  the  a^fi)  satisfy 


the  condition  Q of  Hajek  (1961),  i.e.. 


lim  .,-1/,  . lim/  max  ,,-lr-^N  , . — .21  _ „ 

N-«>  ‘ ^ " 0 ^ N-»~{i<ii<. . .<i  <N  N ^a=l  ^aN -aN  J 0 ' (4,1’ 


V 


9) 


Since,  the  g^  are  real  constants,  by  (4.16)  and  (4.19),  it  follows  by 


some  standard  steps  that 


i£N'V0~l£{lSi<”'i  <N  ■ «• 


20) 


where  d^  = N ^ dj^ (ot)  = 0 . Also,  by  (2.3)-(2.7)  and  (4.16),  it  follows 
that  for  every  g*0,  there  exists  a finite  and  positive  g* , such  that 


N'1I-=lldN(i)-dN]Z  _>  8‘  35  N 


T i 2 


(4.21) 


Thus,  the  d^fi)  satisfy  the  condition  Q of  Hajek  (1961).  Further,  by 


(4-17),  I^jf?  = Ii=1Im=l(cmi)2  = ^i  = l~i%  = p (scc  (4'5))  and  by  (4,2)' 
(4.4),  (4.5)  and  (4.17), 


max  f 2 


l<i<N  i 


-+  0 as  N 00 


(4.22) 


so  that  the  f^  satisfy  the  Noether  condition.  Hence,  the  asymptotic 
normality  of  (4.18)  follows  directly  by  an  appeal  to  the  Hajek  (1961) 


theorem,  and  the  convergence  of  f.d.d.'s  of  (W^  to  W holds. 


Chattcrjec  and  Sen  (1973)  have  established  the  tightness  of  K for 

'iM  t r 

the  special  case  of  p=l.  As  such,  using  their  result  coordinate  wise, 


-14- 


we  obtain  that  for  every  c>0  and  n>0,  there  exists  a 6:  0 < <5  < 1 
and  an  integer  N^,  such  that  for  N > and  each  j(=l,...,p), 

P{sup[|l\'p^  (t)-N^  (s)  | : 0 s s < t s s + 6 < 1]  > ep  2}  < p . (4.23) 

, r h , r 

On  the  other  hand,  | | r(t) -WNjl.(s)  I I 2 = Ijal  [W^^(t) -W^(s)  ]“,  so  that 
by  (4.23), 

P{sup[| |WN  r(t)-WN  r(s)  I I : 0 < s < t s S+6 < 1]  > e}  <n,  V N < NQ  , (4.24) 

while  W (0)  = 0 with  probabi lity  ' 1 (by  definition).  Hence,  the  proof 

'-IN  , J 

of  the  tightness  of  {\V  ) is  complete.  Q.F;..D. 

Let  now  8^  R be  the  signma-field  genci-ated  by  / - (S^,  , . . . ,S^R) , 
OSkSN,  8^  q being  the  trivial  sigma- field.  Then,  for  every  N(>  1), 

8^.  R is  non-decreasing  in  k.  Further,  we  define  as  in  (2.10). 

Lemma  4,2.  Under  HQ  in  (2.2),  {T^  R,  6N  R ; 0 < k < N)  is  a martingale 

for  every  N(>  1),  (i^  R,  B^k;  0<k<N>  'I'S  a submartingale. 

Proof.  Under  H^  in  (2.2),  by  Lemma  4.1  of  Chatterjec  and  Sen  (1973),  it 
follows  directly  that  for  each  j (=l,...,p), 

E(TNjk+llSN  k}  = TNjk  a-C-  for  every  k:  0<k<N-l  , (4.25) 

and  hence,  E (T^  k+jlsN  R)  = XN  R a.e.,  V 0sk<N-l.  The  same  is  true 
for  T,.  . . Further  bv  (2.13),  (2.15),  (2.16)  and  (2.19).  LM  (being 
the  Euclidean  norm)  is  convex  in  R,  and  hence,  the  martingale  pro- 
perty of  {T^  R)  along  with  the  conditional  form  of  the  Jensen  inequa- 
lity yields  the  submart  ingale  property  of  R) . Q.E.P. 
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For  every  t:  0 < t < 1 , let  us  define  now 

f*  ■»  _ i f f ^ ^ 

v(t)  = J ^(u)du  + (1-t)  1 1 j 4>(u)duj  , 

so  that  v(t)  is  in  t c (0,1),  v(0)  = v(t)  = 0 and  v(l)  = 

liin  . . , . 2 

til  1 • 

Lemm3  4.5.  Undci ’ (2.3)-(2.7),  k/N  •>  t : t£  [0,1]  insures  that 


(4.26) 


S*,k  = ^'^i^N.k^  - V(t)  * as  N 


(4.27) 


where  the  ^(i)  are  defined  in  (2.12). 

The  proof  follows  along  the  lines  of  the  proof  of  Lemma  4.2  and 
Theorem  4.2  of  Chatterjee  and  Sen  (1973),  and  hence,  is  omitted. 

Let  us  now  introduce  the  following: 


- V(t)dt  = lp _J\r (t)dt  = lp  (0  , 

Jq  ■>'  Jq  j j 1 j 


(4.28) 


(4.29) 


1 2 

where  the  ok  (=jpWj (t)dt)  are  i.i.d.  nonnegative  r.v.  From  Theorem 
4.1  we  conclude  that  under  the  hypothesis  of  Theorem  4.1, 


M*.  ^ U°  and  K*  % w*  as  N . 

N,r  p TJ,r  p 


(4.30) 


The  characteristic  function  (c.f.)  g(Q)  of  u)^  (or  any  uk  , j > 1) 
is  given  by  [viz.,  Dugue  (1969)] 


g(0)  = (Cos/2i0)  2 = 7 f{l-2i0u.}  u 

k=l  K 


’2;  u.  = 1 j , k>  1 . (4.31) 

K IT  (2k-  1 ) 


I 


L 
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Therefore,  the  c.f.  gP(0)  of  Wp  i-s  given  by 


g°(0)  = [g(0)]P  = THl-SieuJ 
p k=l  K 


■v/: 


(4.32) 


Note  that  if  {if,  j > 1}  be  a squencc  of  i.i.d.r.v.  where  l).  has  the 


central  chi-square  distribution  with  p degrees  of  freedom,  then  for 
every  c. 


E(exp[itcU ] } = [l-2itc]  p/^ 


(4.33) 


so  that  from  (4.32)  and  (4.33),  we  have 

0 V r°°  r a <_2 , , 2i 


to 


= Xk=  i {4/ir^  (2k - i ) ^ }Uk  , 


(4.34) 


where  = stands  for  the  equality  of  distributions.  From  (4.30)  and 


(4.34),  we  conclude  that  under  11  , as  N -*■  00  , 


*£>r  " Ik  = It4/*2(2k-l)2}Uk  = u[J  , say 


(4.33) 


Since,  we  do  not  know  the  distribution  of  Uq  in  any  closed  form,  we 


have  obtained  (by  simulation  studies)  the  empirical  percentile  points: 
these  will  be  reported  in  Section  7. 

For  p=l,  the  distribution  of  (o*  = | (t)  | is  well-known  and 

is  given  by  [viz.,  Billingsley  (p.7 9;  1968)] 


P(u)*  < x}  = ^_ro(-I)k(<5>((2k+l)x)  - <t>((2k-l)x) ] , x>0  , (4.36) 


where  4>  is  the  standard  normal  d.f.  Parallel  expressions  for  p>  1 are 
not  known  and  remain  as  challenging  problems  for  probabi 1 ists.  In  Section 
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7,  we  have  also  dcrivcJ,  through  simulation  studies,  some  empirical 
values  for  the  percentile  points  of  the  d.f.  of  w*  for  p-  4 . 

We  may  however,  note  that,  by  definition,  (Y(t),  tel)  is  a sub- 
martingale, so  that  for  every  x > 0, 

p{£j  v(t)  > x} . r y™  > y}  (e > o) 


inf  / 
= 0>0  \ 


-Ox 

c (1-29) 


-P/2| 


as  Y2(l)  = has  the  central  chi-square  d.f.  with  p degrees 

of  freedom;  note  that  the  inequality  in  (4.37)  is  based  on  the  Kolmogorov 

2 

inequality  for  submart ingalcs.  For  x > p,  the  right  hand  side  reduces 


(4 . 37) 


(e/p)p/2xV’5X  (*SH{/p7 1 2p/2\fT2V\Pe''* 


(4.38) 


On  the  other  hand,  as  oj*  > u>* , we  have  by  (4.36),  for  every  x > 0, 

P{w*  > x}  ? (oj*  > x}  = 1 - P{<  s x) 

P 1 1 

= 2 [ 1 -<t>(x)  ] + 2^”=1  (-l)k_1  [4>(  (2k+l)x) -4>((2k-l)  x)  ] (4.39) 


(~  4 [ 1 -<J>(x)  ] when  x is  not  small) 


From  (4. 37) - (4 . 39) , we  obtain  that 


x-M.rof-2x  "log  P{u>*  > x)  J = 1 , for  every  p ? 1 . 


(4.40) 
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Wc  shall  find  this  result  quite  useful  in  Section  6.  In  passing,  we  may 
remark  that  in  (4 . 37) - (4 . 38) , instead  of  using  the  Kolmogorov  inequality, 
we  could  have  used  the  Doob  uperossing  inequality  for  submartingales  and 
obtain  the  upper  bound 


-Ox  1 


0>°l  2P/2J-/2  J 2 

11  u>-  q 1 og  2+x 

' ^{e-°,‘22(l-2C)-f/22-P/2(t?72)-1 


j' ■ 1c-<.(l-2«u(lu} 


p/ 2—1  -Ijv  1 
v^'  c dv 


v>(l-20)[x - | log  2] 


For  specific  values  of  p (viz.,  p = 2) , the  right  hand  side  of  (4.41) 
can  be  worked  out  explicitly  and  the  same  is  somewhat  sharper  than  (4.38). 


5.  ASYMPTOTIC  NON-NULL  DISTRIBUTION  THEORY 

With  the  intention  of  studying  the  (asymptotic)  power  properties 
of  the  proposed  PCS  tests,  we  proceed  now  to  consider  the  asymptotic  non- 
null distribution  of  K*  ^ and  For  fixed  alternative  hypothesis, 

these  distributions  do  not  exist  and  we  are  left  with  the  task  of  study- 
ing the  rates  of  convergence  of  the  powers  to  1 — as  will  be  done  in  the 
next  section.  On  the  other  hand,  as  is  the  usual  fashion,  we  may  consider 
a sequence  of  local  aternative  hypotheses,  chosen  so  carefully  that  under 
such  a case,  the  asymptotic  non -null  distributions  are  properly  defined 
and  the  powers  are  bounded  away  from  1 . 

Towards  the  end  of  Section  3,  we  have  observed  that  K*  ^ and  r 
are  both  invariant  under  non-singular  transformations  on  the  regression 


t 

1 
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vectors  — hence,  we  nay,  without  any  loss  of  generality,  use  a canonical 
reduction.  We  assume  that  there  is  a triangular  array  (x^,  l<i<N, 

N > 1 } of  rowwise  independent  r.v.,  where  for  each  N, 


< x}  = F(x-30-3'c*i)  , l<i<N  , - 00  < x < 00  , 


(5.1) 


where  c*  = (c* . , . . . ,c*  . ) ' , i = 1, . . . ,N,  and 
~Ni  Nil  Npi  9 > > > 

. c* . = 0 and  TN  .c*  c*'  = I , (5.2) 

‘•i  = l~Ni  ~ *-i=l'~Ni~Ni  ~p  v ’ 

and  in  this  case,  (4.2)  reduces  to  l<j<p|i<i<\'  cN'ji}  As  bcf°rc> 

0=0  and  we  frame  a sequence  Ol^}  of  alternative  hypotheses  by 
letting 

IIN:  (5. 1)-  (5.2)  hold  with  0 * 0 . (5.3) 

Regarding  F in  (5.1),  we  assume  that  it  has  an  absolutely  continuous 
density  function  f with  a finite  Fisher  information 
f°° 

1(f)  = J [f ' (x) / f (x) ] 2dF (x)  , where  f ’ (x)  = ^ f(x)  . (5.4) 


In  addition,  we  define 


<J>°(u)  = -f’(F‘1(u))/f(F'1(x))  , 0 < u < 1 , 


(5.5) 


so  that  /J<j>°(u)du  = o and  [<f>U(u)  ] “du  = 1(f)  < °°.  Also,  for  sc  (0,1), 
we  let 


rlr.O, 


4>s(u)  = 


<J>  (u)  , 0 < u ^ s , 

fl 


(5.6) 


<J>*  = (1-s) 


-1 


4>  ( t ) d t , s < u < 1 , 
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and  define  <}>  (u) , 0<u<l  in  the  same  way.  Let  then  for  se  (0.1), 


p(s)  = 


n ) 1- 

<t>  (u) <p  (u) du  /1(f)  = 


$®(u)<^(u)du 


/IJ(f)  , (5.6) 


and  finally,  we  define  v(s),  0<ssl  as  in  (4.26).  Then,  the  follow- 
theorem  provides  the  bases  for  subsequent  results  of  this  section. 


Theorem  5.1.  Under  (5.1)-(5.4)  and  the  conditions  on  the  score  function 
<p  assumed  in  Section  2,  W , defined  by  (4.7)  [but  for  the  triangu- 
lar  array  of  r.v.  's  in  (5.1)],  converges  in  law  in  the  J ^-topology  on 
D[0,1]  to  a p-variate  Gaussian  function  W + y,  where  W is  defined 
before  (4.9)  and  £ = (ja(t.)  , tel}  is  given  by 

y(t)  = gp(v  I(tv(6))[I(f)/v(6)]?-  , 0<t<l  . (5.8) 

Proof.  Let  P^  and  P^  be  respectively  the  joint  d.f.  of  (X^j , . . . ,X^^) 

under  and  11^  and  let  pj^  and  be  the  same  for  (Z^ Z^) , 

for  k<N.  Then,  by  the  results  of  Chapter  VI  of  Hajck  and  ^idak  (1967), 
we  conclude  that  under  (5.1)-(5.4),  (P  } is  contiguous  to  { p5? } , and 
this  insures  that  (P^.  h<N}  is  also  contiguous  to  {pj^>  k<N}. 

Hence,  we  may  proceed  along  the  lines  of  the  proof  of  Theorem  2 of  Sen 
(1976a)  and  show  that  the  tightness  of  under  Hg,  insures  the 

same  under  the  contiguous  alternatives  •( } . The  convergence  of  f.d.d.'s 
of  K - y to  those  of  W also  follows  by  an  appeal  to  contiguity 
(when  (ll  } holds)  and  the  earlier  part  of  the  proof  of  Theorem  4.1 
along  the  same  line  as  in  Theorem  2 of  Sen  (1976a).  Q.E.D. 
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Recalling  that  W(t)  = (W^(t),...,W  (t))',  tel  and  p(t)  = 


(JJj  (t)  , . . . ,p  }(t))  ' , tel,  we  obtain  from  (>.5),  (3.6),  (4.7),  (4.8) 
and  Theorem  5.1  that  under  the  hypothesis  of  Theorem  5.1,  (under  {ll^}) , 


* 2 sup/  [7-p  . 

S,-,r  tcI\Lfj*l(Wj 


(t)nyt))- 


V 


(5.9) 


V 

-V 


(W. (t)+p. (t) }2dt  . 
0 J J 


(5.10) 


Note  that  the  Pj(t),  defined  hy  (5.8),  are  not,  in  general,  linear  func- 
tions of  t,  and  as  in  Section  4,  the  exact  distributions  for  the  right 
hand  sides  of  (5.9)  and  (5.10)  arc  difficult  to  obtain.  As  a result,  it 
is  difficult  to  draw  more  indepth  conclusions  about  the  relative  perfor- 
mance of  these  PCS  tests  for  contiguous  alternatives.  For  this  reason,  in 
the  next  section,  we  take  recourse  to  the  Bahadur  efficiency,  where  under 
(fixed  but)  close  alternatives,  we  have  some  meaningful  comparisons  of 
the  different  test  statistics  and  score  functions. 


6.  BAHADUR  A.R.E.  OF  PCS  TESTS 

First,  parallel  to  (4.40),  we  derive  a limiting  result  for  the 
tail  probability  of  the  Cramer- von  Mises  type  statistics.  Note  that 
by  (4.35), 

00 

“2  U1  < l,0  = 'TU1  + Ik=2{4/7T2(2k'1)2}Uk  ’ (C’-n 

n TT 

where  each  IK  , j-1,  has  the  central  chi-square  distribution  with  p 
degrees  of  freedom,  so  that  for  every  X>0, 


' 


- i.i~- 


where  c(e)  (<  °°)  depends  on  t(0  < c < ir"/S) . Thus,  choosing  e(>  0) 
arbitrarily  small,  we  obtain  from  (6.5)  that 

2 


lim  inf  P 2 ni-,.0  „ ,27 

X—  L‘?  R 0 _ 


> - 2e  , V e(0  < e < tt2/8)  , (6.6) 


so  that  from  (6.4)  and  (6.6),  we  conclude  that 


limP  2 nr,,0  ^ > 2 P]  2 

1 log  PUI0  > X ) = IT 


X-p  ,2 


/4 


(6.7) 


In  some  specific  cases,  the  stochastic  convergence  of  the  L , , HN 

1\  j K 

can  be  studied  for  general  alternative  [viz.,  Chatterjee  and  Sen  (1973)  for 
the  two  sample  problem].  But,  in  general,  it  demands  extra  regularity  con- 
ditions. Let  us  assume  that.  X^,...,X^  are  independent  with  d.f. 's 
Fj,...,Fn  and  defining  c^  as  in  (4.5),  we  let 

F(N)(x)=N‘1^=iFi(x),  D(n)(x)=N-^=1c.F.(x),  - ~<x<~  , (6.8) 


and  defining  be  (2.15),  we  assume  that 


(i)  3.P"  N *C=  C and  c,,  = c both  exist  , 

v ' N-*®  ~N  ~ N-*»  ~N  ~ ’ 

(ii)  ^ ^(\’)(x)  = F(x)  exists  for  all  x (a.e.)  , 

(iii)  D(iNJ)  (x)  = D(x)  exists  for  all  x (a.e.)  . 


(6.9) 

(6.10) 
(6.11) 


Let  us  also  define  V(t)  as  in  (4.26)  and  let 
rx 


T(X)  = 


h 


<{»(F(y))dn(y)  - (l-F(x)  } 0(x) 


<»(F(x))dF(x))  , 


(6.12) 


C(x)  = (T(x))’(t(x))  , 


- °°  < x < °°  , 


(6.13) 
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i*M  = J^xtCfy)]5  , — » < x < « , 


(6.14) 


C°(x)  = ff  S(y)dv(F(y))]/v(F(x))  , 

U j 


00  < < oo 


(6.15) 


Then,  along  the  lines  of  Section  5 of  Chatterjee  and  Sen  (1973,  p.  41), 
it  can  be  shown  by  some  standard  steps  that  under  (2.3)-(2.6),  (2.21)  and 
(6.9) - (6. 12) , as  N -*■  » , 


Kj,r  ->  C*(F  1(6))/v,5(6)  a.  s . , + 5°(F  1 («))  a.s.  , (6.16) 


LN,r/AN,r  " C(F  1(6))/V(6)  a-s‘ 


(6.17) 


In  particular,  for  the  model  (2.1),  if  the  c^  have  all  .bounded  elements 


and  B is  close  to  0,  then  (6. 16) - (6. 17)  simplify  to 


^N,r  ~~~> 


<J>n(u)du  /v(5) 


| +o(l  !gl !) 


= I(£,£‘1«Kf)/v(fi)]'i[osS6p2ct)}  + 0(!!£ll)  , (6. is) 


where  p(s)  is  defined  by  (5.7)  and  ||B||  = B'B; 


-1  PfF  2 - 

> (B'C  J3)  p“ (F ( 


P (F(y)dv(F(y))  [ I (f)/v  (6) ] + o(| |B| |) 


^‘N,r  (£'££) 


= [(6,C',B)I(f)/v(6)]  p“(u)dv(u)  /v(5)  ♦ o(||B||) 
~ ~ ~ J0 


6.19) 


lN,r  /AN,r  ~~>  (B'C“,3)p2(5)  I (f)/v(5)  + o(||3||) 


(6.20) 


2 2 V 

Finally,  note  that  under  in  (2.2),  by  Theorem  4.1,  y -►  U j 

and  hence  (6.3)  applies. 
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By  virtue  of  Theorem  4.1,  (4.30),  (4.40),  (6.2)  (for  U^) 

and  (6. 16)- (6. 17) , we  are  in  a position  to  adapt  the  Bahadur  efficiency 

results  (viz.,  Puri  and  Sen  (1971,  pp.  122-123),  The  BARE  (Bahadur  ARE) 

of  the  Komogorov-Smi rnov  type  test  relative  to  the  terminal  test  based 

on  L.,  is  given  by 
N , r 
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Let  us  now  confine  ourselves  to  local  alternatives  for  which  (6.18)- 
(6.20)  hold  and  in  this  case,  the  limiting  BARE  reduces  to 

2 ([6  2 9 

e*>T  = T y P (u)dv(u) )/[v(o)p  (6)]  ; (6.: 


’ IT  ' W0 

"IT  * toS6p2<t>]''p2<5)  » • 


u)  J ; 


In  the  context  of  optimality  of  score  functions  for  PCS  rank  tests 
for  simple  regression,  Chatterjee  and  Sen  (1973)  and  Sen  (1976b)  have 
studied  the  optimality  of  <J>^(u),  defined  by  (5.5).  It  follows  from 
(5.4)-(5.7)  that  for  (J)  = (unto  a scalar  constant). 


v(t)  = p(t) 


Uo  1 


9 1 f f ) 9 

- i.  / i i i- 


(u)]*duj/j  [$  (u)  ] “du  , 0<t<l  , 


so  that  we  have  from  (6. 24) - (6. 27) , 


e*  _ = Tt  / 1 2 = 0.S225 

m,  r 

c*  w = 12/it2  - 1.2159 
K,M 

* — i 

eK,T  ‘ 

In  this  case,  we  are  naturally  inclinced  towards  using  the  Komogorov- 
Smirnov  type  tests  on  the  ground  of  the  limiting  BARE.  However,  the 
picture  can  be  different  when  41  * ^ • For  example,  suppose  one  uses 
the  exponential  score  <{'(u)  = -l-log(l-u),  0<u<l,  while  the  under- 


lying distribution  is  logistic.  In  this  case,  c*  reduces  to  15/2tt  = 
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0.7599,  so  that  the  Cramer-von  Mises  type  test  appears  to  have  an 

edge  over  the  Kolmogorov-Smirnov  type  tost.  For  this  example,  e*  = 

M.T 

2 

2 rr  /IS  = 1.5159,  so  that  the  terminal  test  is  also  not  as  efficient  as 

the  Cramer-von  Mises  type  test.  An  opposite  picture  holds  when  one  uses 

the  IVilcoxon  scores  (viz.,  ^(u)  = /IT  (u-1/2),  OSuSl)  while  the  under- 

2 2 

lying  d.f.  is  exponential  - here  e*  = 2n  /35  = 0.5640  and  e*  = 3S/2it  = 

M, l 

1.2665.  These  examples  suggest  that  whereas  the  BARli  e*  j.  or  e£  ^ may 
fluctuate  quite  a bit  for  different  score  functions  and  underlying  d.f.'s, 
c*  T — 1 remains  true  under  quite  general  conditions,  tending  to  advocate 
the  use  of  Komogorov-Smirnov  type  of  PCS  tests. 

7.  SIMlILATF.il  PERCENTILES  OF 

NULL  DISTRIBUTIONS  MA  AND  K* 

N,r  N,r 

The  distributions  of  M*  f and  r have  been  shown  to  converge 
weakly  to  some  functionals  of  the  standard  Wiener  process  under  the  null 
hypothesis  and  to  those  of  the  drifted  Wiener  process  under  contiguous 
alternatives,  under  certain  regularity  conditions.  As  we  have  mentioned 
in  Section  4,  the  null  distributions  of  these  processes  are  not  available 
in  workable  form.  We,  therefore,  derive  in  this  section  a few  percentile 
values  of  these  distributions  empirically  through  simulation  studies. 

Consider  n independent  observations,  Y from  the  stan- 

dard normal  distribution.  Let 

k 

sk  = l Yk  , isksn  , 
i=l 


Sq  = 0 by  convention  . 


(7.1) 
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We  define  the  stochastic  process  W^  = (|V  (t),  t c [0,1]}  by  letting 


-"'V) 


(7.2) 


where  n(t)  = max{k:  k^tn}. 

We  note  that  the  sample  process  W is  right  continuous  with  the  left- 
hand  limit  and  hence  belongs  to  the  metric  space  D[0,1]  with  the  pro- 
perties 

EWJt)  = 0 , 

EW*(t)  = n"1 [nt]  , (7.5) 

EKn(t)1Vn(t,)  =n'1[n(tAt,)]'  t't'e  ’ 

wliere  [nt]  and  [n(tAt')]  denote  integral  parts  of  nt  and  n(tAt') 
respectively.  The  maximum  jump  of  the  process  is  given  by 

»,x  i\i 

1SkSn  xS  Ji. 


0 a. s. , as  n 


(7.4) 


Consequently,  as  n gets  large  the  process  W has  a continuous  sample 
path  a.s.  and  has  the  structure  of  the  standard  Brownian  motion  process 
W = (W(t),  t e [0,1]}. 

Let  now  [{Y„K_j,  i = l,...,p]  be  independent  random  samples  each 
of  size  n from  W(0,1).  Then  by  Donskcr's  theorem  [Billingsley  (pp. 
6S-77 ; 1963)] 
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max 

0<krn 


L 7 

l w:n(k/ni 


i-1 


U)* 


1 p n'1  ? 
n I I W.n(k/n) 


i=l  k=0 


V 0 

-*■  U) 


(7.5) 


n * 


wh 


0 


ere  W.  (k/n)  = n 1 Y Y.  . , k = 0,. . . ,n  and  i = 1 , . . 
in  >,11  ’ 

j = l 


,p  and  to*  and 
* P 


are  given  by  (4.28)  and  (4.29)  respectively. 


For  purposes  of  the  simulation  studies,  we  have  generated  the  stan- 
dard normal  deviates  by  using  IBM  scientific  subroutine  GAUSS.  The  sam- 
ple size  n has  been  taken  to  be  1000  and  the  empirical  1 (null)  distri- 
butions of  the  two  processes  have  been  derived  through  1000  independent 
repetitions.  For  details,  see  Majumdar  (1976).  In  Table  .1  and  Table  2 
below,  we  have  furnished  a few  simulated  values  of  the  right  tails  of 
the  two  distributions.  For  p=l,  we  have  given  exact  percentile  values 


of  to*  by  using  the  approximation 


> xj  = 4(1  -4>(x) ) , 


(7.6) 


where  4>(x)  is  the  probability  integral  of  the  standard  normal  distribu- 
tion. 


TAR1.F.  1:  Simulated  values  of  the  null  distribution  of  to  for 

P 


selected  values  of  p and  a 


1* 

'J 

C 

3 • 

4 

.01 

2.81 

3.22 

3.71 

3.89 

.05 

2.24 

2.70 

3.05 

3.31 

.10 

: t -=.x  --  r-r  - 

1 .96 

is 

2. -s 

3 . 04 

--  r s -r- 

*F.xact 
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TABLH  2:  Simulated  values  of  the  null  distribution  of  for 

selected  values  of  p and  a P 


a ~~P  ' 

1 

2 

3 

4 

.01 

2.87 

4.06 

5.57 

6.54 

.05 

1.66 

2.67 

3.49 

4.33 

.10 

1.19 

2.01 

2.80 

3.64 

8.  SOME  GENERAL  REMARKS 


As  has  been  mentioned  in  Section  1,  our  model  (2.1)  includes  the 

multi-sample  location  model  as  a special  case'.  Let  X..^,  1 ^ i ^ n . be 

i.i.d.r.v.  with  a continuous  d.f.  Fj (x) , for  0 s j £ p(^l)  and  let 

N = n_  + . . . + n . Rewriting  X..  = X. , U i i n„  and  X.  . = X , , . , 

0 p ° Oi  l 0 1J  11-+...+11.  ,+i 

^ J 0 j-1 

1 £ i £ iu  , 1 £ j £ p,  and  assuming  the  conventional  location  model  where 
F_.  (x)  = F(x-0j ) , 0 < j < p,  we  observe  that  (2. 1)  holds  with  3^.  = G_.  - Gq. 
l£i£p,0  =0  and  c = . . . = c = 0,  c , = ...=  c = (1 ,0, . . . ,0) ' , 

J 1 ’ 0 0 ~1  ~nQ  ~ ~i1q+1  ~nQ+ni 

. . . ,c  , = ...=  cM  - (0, . . . ,0, 1) ' . The  null  hypothesis  H in 

~n,+...+n  . + 1 ~N  0 

0 p-1 

(2.2)  insures  that  F^ = ...  = F . If,  we  assume  that  the  sample  sizes 


n.,...,n  satisfy  the  conditions 
O p 


A*™  N-1n . = A. : 0 < X.  < 1 , V 0<  j <p  , 
N-*»  j j J J 1 


then  £n  -*•  (Xg,...,X  )'  as  N -*■  00  and,  by  (2. IS),  N !£N  ■*  ‘ 

A.X  ))  , so  that  (2.16)  holds.  Thus,  the  proposed  PCS  rank 

K q K , (J=  1 , . . . , p 

tests  apply  to  the  multi-sample  location  problem  as  well. 
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ln  order  to  test  homogeneity  of  k(=p+l)  samples  for  right-censored 
data  (fixed-plan  censoring),  with  the  smallest  r out  of  N observations 
of  the  combined  sample  being  considered,  Basu  (1967)  has  studied  a genera- 
lized version  of  the  Kruskal -Wallis  test.  The  asymptotic  chi-square  dis- 

2 2 

tribution  of  his  statistic  (cqivalent  to  our  L.,  /A..  for  his  particu- 

1 N , r N , r 

lar  model)  follows  readily  from  our  Theorems  4.1  and  5.1.  In  his  scheme, 

early  termination  of  experimentation  (prior  to  the  r-th  order  statistic) 

has  not  been  advocated,  while  in  our  PCS  procedures,  this  is  no  problem. 

One  can  use  K*  or  M*  . The  BARI-  results  of  Section  6 suggests  that 
N,r  N,r  bb 

2 2 

using  K*.  instead  of  L„  /A.,  allows  an  earlv  termination  without 
N , r N , r N , r 1 

any  loss  of  the  asymptotic  efficiency. 

Suppose  now  that  instead  of  a preassigned  number  r of  failures, 
the  experiment  is  designed  to  continue  at  most  for  a period  of  y time- 
units.  Then,  r(v),  the  number  of  failures  occurring  in  the  time-period 
y,  is  itself  a (non-negative  integer-valued)  random  variable.  As  in 
Section  2 of  Chatterjee  and  Sen  (1973),  the  distribution  theory  of 
Mj\|  r(y)  or  r(y)  Can  developed  [under  the  null  hypothesis  11^ 
in  (2.2)]  under  a conditional  setup,  given  r(y)  =r.  However,  in  prac- 
tice, this  conditional  argument  requires  some  knowledge  on  the  distribu- 
tion of  r(y)  so  that  the  stochastic  limit  of  N *r(y)  is  fairly  known 

in  advance  of  experimentation  (as  the  same  is  needed  to  define 
-) 

A7,  for  both  the  PCS  tests) . We  may  surmount  this  problem  by  work- 
N,r 

ing  with  an  upper  bound  for  N r(y)  [allowing  chance  fluctuation], 
whenever  feasible.  In  the  Department  of  Biostatistics,  University  of 
North  Carolina,  Chapel  Hill,  a seven-year  project  on  the  effect  of  high 
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cholestcrol  on  the  risk  of  heart  attack  is  under  study.  Male  patients 
(over  the  age  35)  are  randomly  allocated  to  either  the  control  or  treat- 
ment groups  and  the  survival  patterns  of  the  two  groups  arc  being  pro- 
gressively studied.  From  independent  sources  (viz.,  U.S.  Life  tables), 
the  seven  year  mortality  rate  for  the  particular  age-pattern  is  roughly 
known  to  be  about  ]\°o,  so  that  for  a sample  of  size  N,  an  upper  (OS* 
or  99%)  confidence  limit  can  be  set  on  the  actual  number  of  failures  in 
this  study  period,  and  with  that  upper  limit,  we  can  set  our  proposed 
PCS  tests.  This  procedure,  though  a bit  conservative,  performs  quite 
well  (in  scope  as  well  as  in  performance)  as  compared  tc  some  parametric 
tests  based  on  particular  forms  of  failure  distributions. 
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